Quarkonium suppression as a probe of a 
saturated gluon plasma? 

V. P. Gongalves * 

Instituto de Fisica e Matemdtica, Univ. Federal de Pelotas 
Caixa Postal 354, 96010-090 Pelotas, RS, BRAZIL 

Abstract: 

A dense parton system is expected to be formed in the early stage of rel- 
ativistic heavy-ion colhsions at RHIC energies and above. The probabihty of 
a quark gluon plasma production and the resulting strength of its signatures 
depends strongly on the initial conditions associated to the distributions of 
partons in the nuclear wave functions. At very high energies, the growth of 
parton distributions should saturate, possibly forming a Color Glass Con- 
densate, which is characterized by a bulk momentum scale Qs- As a direct 
consequence, the possible signatures of the QGP should be Qs-dependent if 
the saturation scenario is valid for RHIC and LHC. In this Letter we assume 
the saturation scenario for the QGP formation and estimate the saturation 
scale dependence of quarkonium suppression. We conclude that, if this sce- 
nario is vahd, the T suppression only occurs at large values of Qs- 
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High energy heavy-ion colhsions offer the opportunity to study the prop- 
erties of the predicted QCD phase transition to a locally deconfined quark- 
gluon plasma (QGP) A dense parton system is expected to be formed 
in the early stage of relativistic heavy-ion collisions at RHIC (Relativistic 
Heavy Ion Collider) energies and above, due to onset of hard and semihard 
parton scatterings. These processes happen in a very short time scale af- 
ter which a dense parton gas will be formed, not immediately in thermal 
and chemical equilibrium. Secondary interactions among the produced par- 
tons may lead to thermalization if the interactions are sufficiently strong, 
allowing to study the subsequent dynamical evolution of the system using 
hydrodynamic-based models 0. Exactly how this parton gas equilibrates 
and a thermalized plasma of quarks and gluons is formed requires solving a 
kinetic equation, given appropriate initial conditions and a detailed knowl- 
edge of the microscopic processes by which the quanta exchange energy and 
momentum. This topic is currently under intense debate p, ^, |^. 

The search for experimental evidence of this transition during the very 
early stage of the AA reactions requires to extract unambiguous character- 
istic signals that survive the complex evolution through the later stages of 
the collision. One of the proposed signatures of the QCD phase transition is 
the suppression of quarkonium production, particularly of the J/V' @. The 
idea of suppression of cc mesons J/ip, etc., is based on the notion that cc 
are produced mainly via primary hard collisions of energetic gluons during 
the preequilibrium stage up to shortly after the plasma formation (before the 
initial temperature drops below the production threshold), and the mesons 
formed from these pairs may subsequently experience deconfinement when 
traversing the region of the plasma. In a QGP, the suppression occurs due 
to the shielding of the cc binding potential by color screening, leading to the 
breakup of the resonance. The cc {J/ip,ip', ...) and bb (T, T', ...) resonances 
have smaller radii than light-quark hadrons and therefore higher tempera- 
tures are needed to dissociate these quarkonium states. 

The probability of a thermalized QGP production and the resulting strength 
of its signatures strongly depends on the initial conditions associated to the 
distributions of partons in the nuclear wave functions. At very high energies, 
the growth of parton distributions should saturate, possibly forming a Color 
Glass Condensate [For a pedagogical presentation see Ref. 0], which is 
characterized by a bulk momentum scale Qg. This regime is characterized 
by the limitation on the maximum phase-space parton density that can be 
reached in the hadron wave function (parton saturation) and very high values 
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of the QCD field strength squared oc l/a^ @]- In this case, the num- 
ber of gluons per unit phase space volume practically saturates and at large 
densities grows only very slowly (logarithmically) as a function of the energy 



10| . If the saturation scale is larger than the QCD scale Aqcd, then this 
system can be studied using weak coupling methods. The magnitude of Qs 
is associated to the behavior of the gluon distribution at high energies, and 
some estimates has been obtained. In general, the predictions are ~ 1 
GeV at RHIC and ~ 2 - 3 GeV at LHC 



Recently, Mueller has proposed ideahzed initial conditions for a satu- 
rated nuclear wave function and shown that the evolution of single particle 



distributions could be described by a nonlinear Landau equation |T3[. This 
equation was numerically solved in Ref. ^ and in a "self-consistent" re- 
laxation time approximation in Ref. 0. One of the main conclusions of 
these studies is that the equilibration time and the initial temperature of the 
plasma have a strong functional dependence on the initial gluon saturation 
scale Qg. As a direct consequence, the possible signatures of this saturated 
gluon plasma should be Q^-dependent if the saturation scenario is valid for 
RHIC and LHC. We denote saturated gluon plasma, the plasma formed from 
the collision of two nuclei with saturated nuclear wave functions. Our goal 
in this Letter is assuming the saturation scenario for the QGP formation, 
estimate the quarkonium suppression and its dependence in the saturation 
scale at LHC energies. 

Lets start from a brief review of the quarkonium suppression in heavy-ion 
collisions. The heavy quark pair leading to the quarkonia are produced in 
such collisions on a very short time-scale ~ l/2mQ, where mq is the mass 
of the heavy quark. The pair develops into the physical resonance over a 
formation time Tf . For a confined surrounding, the bound states of cc and 
bb interacting via two forces: a linear confining potential and a color-Coulomb 
interaction. In the plasma phase, the linear potential is absent due to the 
high temperature leading to deconfinement and the color charge is Debye 
screened by a cloud of surrounding quark-antiquarks pairs which weaken the 
binding force between the QQ pair, thus reducing the color charge seen by 
the other (anti)quark. Basically, above the critical temperature Tc, we have 

V{r) = \- ar ^ V{r) = exp(— /i^ir) 

where a and a (the string tension) are phenomenological parameters and 



is the Debye mass. As discussed in Ref. [14], for values of the screening mass 
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above a certain critical value /xf** {i = J/ilj,ip':Xc,'^,'^',Xb), the screening 
become strong enough for binding to be impossible and the resonance no 
longer forms in the plasma. As a result, the Q and Q drift away from 
each other, after leaving the interacting system to form a charm and bottom 
mesons, which subsequently decays into leptons to be detected. The critical 
screening mass for quarkonia are summarized in Table 1. 

The discussion above implies a procedure to determine if quarkonium 
suppression is expected. Basically, if we known the evolution of the screening 
mass during the collision process we can estimate what quarkonia will be 
effectively suppressed. At screening mass values greater than fif^'^ for values 
of time above (or equal) to the formation time of the bound state (tf), the 
bound state does not be formed. Therefore, we can analyze the quarkonium 



suppression considering the dynamical evolution of the screening mass |T5|. 
An expression currently found in the literature for the screening mass is the 
following: 

/.2 = -^_lim fd'k^lf.'^^fit) , (1) 

where is the strong coupling constant and f{k) is the phase space density 
of gluons. This expression was derived in Ref. |]TB[, following the standard 
calculation (in the Coulomb gauge) of screening in the time-like gluon prop- 
agator in a medium of gluon excitations. In the case of an ideal equilibrium 
gluon gas, f{k) is the Bose- Einstein distribution and the above expression 
implies the well-known screening mass /i^ = ATragT"^- Although the extrap- 
olation from this expression for a system of highly non-equilibrated gluons 
produced in the early stage of heavy ion collisions be useful in the literature, 
is but may not be well defined. 

An alternative expression for the screening mass has been derived in Ref. 



T3| considering the saturation scenario for the QGP. This scenario is based 
on the observation that at very high energies the parton densities in large 
nuclei reach saturation and the number of partons becomes very large, al- 
lowing a description by a classical effective field theory (EFT) [0. In this 
approach (the McLerran-Venugopalan approach), the valence quarks in the 
nuclei are treated as classical sources of the color charges, with the gluon 
distribution for very high nuclei obtained by solving the classical Yang-Mills 
equation. The behavior of the gluon distribution is characterized by the sat- 
uration momentum Qg, below which the distribution reach their maximum 
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value, i.e., the distribution exhibits saturation. When quantum corrections 



are included, a Wilsonian renormalization group picture emerges [|T8[. The 
structure of classical fields remains intact, the scale Qs grows because hard 
gluons (harder than the x scale of interest) in the source increases the typical 
size of color fluctuations. Due to the large occupation number, the small x 
partons described by the EFT constitute a color glass condensate 0. Since 
the classical fields of the two nuclei are known, they provide the initial condi- 
tions for the gluon fields produced in the nuclear collision, with the particle 
production after the colhsion described by the space-time evolution of classi- 
cal gauge fields. The problem of finding the evolution of the classical gauge 
fields in nucleus-nucleus collisions has been addressed in lattice simulations 
|]TP[, allowing to compute the energy and the number of produced gluons at 



central rapidities pO|. On the other hand, in Ref. JTSf Mueller has investi- 



gated the approach to equilibrium in the McLerran-Venugopalan approach, 
considering only elastic scattering and derived a Landau-type equation for 
the single particle distributions f{x,p). The idealized initial condition for 
this transport equation is assumed as given by 

fix,p) = ^^-s{pMQl-Pt) , (2) 

asJMc t 

where c ~ 1.3 is the parton liberation coefficient accounting for the transfor- 
mation of virtual partons in the initial state to the on-shell partons in the 
final state. In Ref. the transport equation has since been solved nu- 
merically using the above initial condition and the approach to equilibrium 
studied quantitatively, verifying that the equilibration time and the initial 
temperature of the saturated gluon plasma have a strong functional depen- 
dence on the initial gluon saturation scale Qg- Furthermore, Mueller has 
assumed that at early times, in the maximum distance 1 / g^m corresponding 
to the minimal momentum transfer g^m one has at most one scattering. This 
condition regularize the logarithmic coUinear divergence in the transport in- 
tegral and provides an alternative expression for the dynamical screening in 
the system of gluons produced at central rapidities, which are completely out 
of equilibrium. One obtains W^ 



2 ,CasA^c^2 Qs 



- (^r- . (3) 

As shown in Ref. P], this expression for the dynamical screening mass is 
larger than the equilibrium expression Eq. (|lD for early times and is almost 
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identical at late times [The difference between the two expressions is ~ 0.05 
GeV for t ~ 0.04/m and is approximately zero for t = Ifm]. As the heavy 
quark production and the formation of bound states are predicted to occur 
in the initial stage of the deconfined phase, we will use the expression (|^) for 
the dynamical screening mass in our analyzes, which implies an upper limit 
for the predictions of quarkonium suppression. 

Now let us consider a central collision in a nucleus-nucleus collision, which 
results in the formation of a quark-gluon plasma. Consider the central slice 
{y = 0) of a collision which leads to the production of a heavy quark pair 
from initial fusion and to a QGP. During the initial interaction period, heavy 
quarks will be produced through gluon fusion and quark- antiquark annihila- 
tion. Like gluon and light quark production, heavy quark production through 
the initial fusion is very sensitive to the parton distributions inside the nuclei. 
Although the cross sections for heavy quark production in RHIC and LHC 



are strongly modified by the presence of the high density effects [21, 22|, the 



quarkonium production rate should be large enough for any suppression to 
be observable. The charm quarks will be produced on the order of a time 
scale l/2mc ~ 0.07/m, while the bottom quarks will similarly be produced 
over a time of ~ 0.02/m. Thus immediately upon production, these quarks 
will find themselves in a deconfined matter which is highly out of equilibrium 
and rapidly evolving towards kinetic equilibration. In the plasma rest frame. 



the QQ forms a bound state in the time tp = 'jTp, where 7 = y 1 + {pr/My 
[M is the onium mass) and the formation times, Tp, are = 0.76/m and 
Tp^^ = 0.89/m. Here we will concentrate in the small pt region, where 
tp ^ Tp. The color charge of the heavy quark will be subject to screening 
due to the presence of the saturated gluon plasma. If the dynamical screen- 
ing mass fi [Eq. (|^)] is larger than the critical screening mass fif'^'^ [Table 1], 
the Q and Q cannot form a bound state and the QQ system will dissociate 
into a separate Q and Q, which subsequently hadronize by combining with 
light quarks or antiquarks to emerge as "open flavors" mesons. We will con- 
centrate our analyzes in the T suppression at LHC energies, since at current 
energies, the situation of J/ip suppression is rather ambiguous because the 
bound state can also break up through interactions with nucleons and co- 
moving hadrons, i.e., QGP production has not been proved to be the unique 
explanation of the observed suppression even though an increased density 
of secondary production is needed. Because the T is much smaller than cc, 
a much higher dynamical screening mass is necessary to dissociate the T, 
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providing a valuable tool to determine the initial state of the system and the 
characteristics of the plasma p^. Moreover, we will consider the dynami- 
cal evolution of the screening mass for t > to, where to is inferior limit for 
the application of the transport theory. For the LHC energies which we are 
interested to < 0.18/m If. 

In Fig. we present the time dependence of the dynamical screening 
mass as predicted in a saturated gluon plasma for some values of the initial 
gluon saturation scale. We can see that our predictions are very sensitive on 
the saturation scale. Considering that the coupling constant is constant and 
equal to 0.3, we have that for small values of the saturation scale {Qs ~ 1.0 
GeV) our results predict that the meson T is not dissociated, causing it to 
escape unscathed. For medium values of Qs (~ 2.0 GeV) the dynamical 
screening mass is larger than the dissociation mass only for the early stage 
of evolution, where the formation of a bound state is unexpected. A similar 
scenario is obtained if we consider the values of c (=1.23), as (=0.6) and Qs 
(=1.41) derived in Ref. [0], where the recent data of the PHOBOS collabo- 
ration for the charged multiplicity distribution were analyzed considering the 
saturation scenario. Therefore, if the saturation scenario is valid at RHIC 
with the typical values for ag and Qs extract from the multiplicity distribu- 
tion we does not expect the T suppression for the current energies in this 
collider {^/s = ISOGeV). Only for large values of the saturation scale we 
predict the T suppression, since by the time that the upsilon is formed the 
dynamical screening mass is above of /i^*^? required to inhibit its formation. 

In Fig. ^ we present the dependence of the dynamical screening mass in 
the the coupling constant for two typical values of the saturation scale. We 
analyze the behavior of the dynamical screening mass in the interval of times 
between the inferior limit for the application of the transport theory {t ^ 0.2 
fm) and the characteristic value of formation time {t ~ 0.8 fm). We verify 
that if the plasma is characterized for a small value of Qs , the suppression 
of the meson T is only expected for early times of evolution and large val- 
ues of the coupling constant, where the perturbative calculations becomes 
questionable. For large values of the saturation scale, the T suppression is 
predicted for as > 0.3 in the initial stage of the kinetic evolution. However, 
at late stages the T dissociation is not expected. For comparison, we also 
present in Fig. |^ the dissociation mass of the J/\E' meson. We verify that in 
this case the complete suppression is expected to occur in the presence of a 
saturated gluon plasma. 

In Fig. ^ we present the dependence of the dynamical screening mass 
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in the saturation scale for a typical value of the coupling constant. We 
verify that the dissociation of the T must occur only if the saturated gluon 
plasma is characterized by large values of the momentum saturation scale. 
In other words, the presence of T suppression indicates that if a saturation 
scenario is valid, the initial nuclear wave function are characterized for a 
bulk momentum scale Qg > 1.5 GeV. Therefore, the T suppression may be 
an indirect probe of the nuclear wave functions, and consequently, from the 
Color Glass Condensate, allowing to constraint the initial gluon saturation 
scale. 

Some comments are in order. In this Letter we have assumed that the 
saturation scale is in the range 1-3 GeV at the RHIC and LHC energies, and 
derive the respective pattern of T suppression. A more precise constraint of 
the saturation scale value can be obtained if we consider the eA processes, 
which could be analyzed in the future eA colliders at RHIC and HERA. In 
particular, the analyzes of the slope of the nuclear structure function |[T^ and 
the ratio between the transverse and longitudinal structure function |25| are 



potential observables to determinate the value of Qg- Our results demonstrate 
that this independent constraint of the saturation scale is fundamental for 
the description of the observables in AA process. 

Lets compare our results with other predictions of T suppression. In Ref. 
T3| the expression (|I|) has been used to analyze the evolution of the dynam- 



ical screening, as well as the onset of quarkonium suppression, considering 
initial conditions provided by the parton cascade model. One of the main 
conclusions of that analyzes is that the T suppression is only expected at 
LHC energies. Distinctly from our results, that predicts that the suppression 
must occur only for the quarkonia produced in the initial stages and at large 
values of the saturation scale, these authors predicts the complete suppres- 
sion, since the screening mass is ever above of the dissociation mass in the 
dynamical evolution. Such result is associated to the initial conditions used 
as input in the calculations, which does not consider a saturation criterium 
in the nuclear wave functions. Similar situation occur if we compare our 
results with that obtained in Ref. where the quarkonium suppression 



in an equilibrating quark-gluon plasma was evaluated considering two mech- 
anisms: the Debye screening and the gluonic dissociation, concluding that 
the Debye screening is not effective for T suppression at the LHC energy. 
This result derives from the self screened parton cascade model, which leads 
to small Debye mass and implies that at the formation time of the T, the 
Debye mass drops below the value of fi^^^, required to inhibit its formation. 
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causing it to escape unscathed. Therefore, the main difference between the 
predictions for T suppression is associated to the use of distinct initial condi- 
tion for the dynamical evolution. Consequently, independent searches of the 
gluon saturation of the nuclear wave function are required to discriminate 
between the distinct approaches. 

In brief, we have seen that if the saturation scenario is valid for RHIC and 
LHC, with the nuclear wave functions characterized by a bulk momentum 
scale Qs, the pattern of quarkonium suppression is also Qs dependent. We 
have verified that the T suppression is only expected at large values of the 
saturation scale, allow us to constraint the values of this scale. In other 
words, the presence of T suppression is an indicative of the presence of a 
saturated gluon plasma, characterized by a large value of the initial gluon 
saturation scale in the nuclear wave functions. Of course, more quantitative 
analyzes can be made if this scale was derived in an independent process. 
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Tables 



State 






Xc 


T 






lif' (GeV) 


0.70 


0.36 


0.34 


1.57 


0.67 


0.56 



Table 1: The dissociation screening mass for quarkonia. 
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Figure Captions 



Fig. |T]: The time dependence of the screening mass in a saturated gluon 
plasma for some typical values of the saturation scale Qs and the coupling 
constant a^. The horizontal line characterize the dissociation screening mass 
for the meson T. 

Fig. ^: The dependence of the screening mass in the coupling constant for 
two instants in the saturated plasma evolution. The long dashed and dotted- 
dashed curves characterize the dissociation screening mass for the mesons T 
and J/\l/, respectively. 

Fig. 1^: The dependence of the screening mass in the saturation scale 
for two instants in the saturated plasma evolution. The long dashed curve 
characterize the dissociation screening mass for the meson T. 
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